Theory of Spin Exciton in the Kondo Semiconductor YbBi 2 
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The Kondo semiconductor YbBi2 exhibits a spin and charge gap of approximately 15 meV. 
Close to the gap energy narrow dispersive collective excitations were identified by previous inelastic 
neutron scattering experiments. We present a theoretical analysis of these excitations. Starting from 
a periodic Anderson model for crystalline electric field (CEF) split 4f states we derive the hybridized 
quasiparticle bands in slave boson mean-field approximation and calculate the momentum dependent 
dynamical susceptibility in random phase approximation (RPA). We show that a small difference in 
the hybridization of the two CEF (quasi-) quartets leads to the appearance of two dispersive spin 
resonance excitations at the continuum threshold. Their intensity is largest at the antiferromagnetic 
(AF) zone boundary point and they have an upward dispersion which merges with the continuum 
less than halfway into the Brillouin zone. Our theoretical analysis explains the most salient features 
of previously unexplained experiments on the magnetic excitations of YbBi2. 

PACS numbers: 71.27.+a, 75.40.Gb, 71.70.Ch 



The so-called Kondo insulators or semiconductors like, 
e.g., CeNiSn, SmBg and YbP>i2 represent a special class 
of strongly correlated electrons [l|. In these compounds 
the conduction electrons hybridize with nearly localized 
4f electrons. The Coulomb repulsion of the latter results 
in a small energy gapQ of order 10 meV at the Fermi 
leveldi. 

At temperatures higher than the gap energy these ma- 
terials behave like Kondo metals exhibiting their typical 
spin fluctuation spectrum. But a low temperatures a 
spin and charge gap opens indicating the formation of an 
insulating singlet ground state 0, Hf- This may be con- 
cluded from the total suppression of the local moment 
in the susceptibility and from the semiconducting behav- 
ior of the resistivity, respectively [f|. The gap formation 
may also be seen directly in the dynamical susceptibility 
and finite frequency conductivity as probed in inelastic 
neutron scattering (INS) and optical conductivity exper- 
iments. In cubic YbBi2 the spin 0] and charge [1] gap 
obtained in this way are approximately equal to 15 meV 
but in general they need not be identical. 

In addition unpolarized Q and polarized 0] INS has 
found an interesting dispersive fine structure around this 
threshold energy. Three excitation branches have been 
identified with energies 15, 20 and 38 meV, respectively 
by analyzing the spectral function of the dynamical sus- 
ceptibility. Since the lower two INS peaks are narrow 
and mostly centered at the zone boundary L-point with 
Q = (tt, 7t, 7t) they may be associated with the forma- 
tion of a collective heavy quasiparticle spin resonance 
exciton appearing around the spin gap threshold @, @] 
and driven by heavy quasiparticle interactions. The col- 
lective modes remain visible in the 20 meV region up 
to T = 159 K 0, Hlj |. Similar spin resonance phe- 
nomena appear as result of feedback effect in unconven- 



tional heavy-fermion superconductors below the quasi- 
particle continuum threshold at 2Ao where Ao is the 
gap amplitude [l2T|. The upper peak is much broader and 
shows little dispersion. It is also rapidly suppressed with 
increasing temperature. It has been associated with con- 
tinuum excitations [l3} also visible in a broad maximum 
in the optical conductivity [1] around 38 meV. 

Furthermore INS experiments on Ybi_ x Lu x B 12 com- 
pounds for different Lu concentrations have indicated 
that the disruption of coherence on the Yb sublattice 
primarily affects the narrow peak structures occurring 
at 15 — 20 meV in pure YbBi2 compound, whereas the 
spin gap and the broad magnetic signal around 38 meV 
remain almost unchanged 14 1 . 

These intriguing experimental results have commonly 
been interpreted in a qualitative way within the spin ex- 
citon scenario @, H, EH but an alternative model was also 
proposed [l6| . However no analysis of the former was at- 
tempted sofar although it is of fundamental importance 
to understand the microscopic origin and fine structure 
of the spin gap in Kondo semiconductors. In this com- 
munication we show in detail how the spin exciton bands 
in YbBi2 arise on the background of a single-particle 
continuum at the spin gap edge. We discuss the origin of 
the splitting into two modes, its connection to CEF ef- 
fects as well as their spectral shape and dispersion. Our 
investigations clarify the underlying microscopic physics 
of these intriguing and for a long time unexplained ob- 
servations. 

Our starting point is the hybridization-gap picture 
based on the periodic Anderson model which is the most 
widely accepted for the description of Kondo semicon- 
ductors. Using the mean-field slave boson approximation 
for CEF split 4f states of Yb we calculate the hybridized 
bands. With an empirical model for the quasiparticle 
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FIG. 1: Dynamical susceptibility in direction (1,1,1) for q = 

(direct gap) and q = Q (indirect gap) versus energy, for pjG. 2: Contour plot of real part of noninteracting dynamical 
degenerate quasiparticle bands with Vi = It and SV = susceptibility for degenerate bands with Vi = It and 8V = 
(V = 0.3CK, Jq = 0). Inset shows the density of states for ^ = Q ^ Jq = Q) in the direction (x,l,l)- 
two CEF-split quasiparticle bands (e^ = 0.08i and b = 0.41). 
The green curve corresponds to the band V\ — t and the red 
one to the band V 2 = Vi + 0.131 

hole Hilbert space is then 



interactions we evaluate the momentum dependent dy- 
namical magnetic susceptibility in RPA. Its imaginary 
part is proportional to the INS spectrum. We obtain 
sharp resonance features around the continuum thresh- 
old and wave vectors not too far from the zone boundary 
L-point. Away from this point the resonance peaks dis- 
perse upwards in energy and broaden. They merge into 
the single particle continuum less than halfway into the 
Brillouin zone (BZ), which describes the basic experi- 
mental facts. In addition our calculation suggests that 
CEF splitting and associated CEF orbital dependence of 
hybridization are responsible for the observed splitting 
into two dispersive resonance modes. 

The Yb electronic configuration is 4/ 13 correspond- 
ing to a single hole in the 4/-shell 17] . Therefore we 
consider the Anderson lattice model with a /-hole in 
a j = 7/2 state, including the CEF effect, i.e., H t = 
Hf + H d + Hf_ d + H c . Here Hf describes the lattice of 
the localized, CEF-split 4/-holes, H d the conduction elec- 
trons and Hf-d is the hybridization between both. Fi- 
nally He is the Coulomb interaction with an on-site hole 
repulsion f///. Our model assumes the limit [///—> oo 
where doubly occupied (hole) states (4f 12 ) are excluded 
and the two possible Yb configurations are either 4/ 14 
or 4/ 13 . The one without a 4/ hole, i.e., 4.f 14 can be 
accounted for by an auxiliary boson b\ [18]. In cubic 
symmetry the j = 7/2 multiplet is split by the CEF into 
a quartet Tg ground state and two excited doublet states. 
The latter may be treated as a quasi quartet T' s according 
to INS results at higher temperatures [3]. The two quar- 
tets (index T = 1,2) have energies Ai = and A 2 > 0. 
The model Hamiltonian in the restricted zero- and one- 



H = $>/ + A 7 )/; 7 / i7 + ^e k 4 7 d k7 

+ N-^^V^e^fld^h + c.c), (1) 
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Here 7 = (r, m) where r = 1, 2 denotes the quartet and 
m = 1 — 4 is the orbital degeneracy index. Furthermore 
the local constraint Qi = b\bi + ^ /2,/i 7 = 1 nas to 
be respected for all i. Therefore the total Hamiltonian 
including the constraint is H — Ab Y^iiQi ~ 1)' wnere 
is the Lagrange multiplier. Here the fj create /-holes at 

lattice site i in CEF state 7, and the create the holes 
in the conduction band with wave vector k and CEF state 
index 7. The f-orbital energy is while A 7 — Ar is the 
CEF excitation energy, and N s is the number of lattice 
sites. Finally is the hybridization energy between 4f 
and conduction holes. In the following, the k dependence 
of the hybridization energy is neglected, i.e., Vk 7 = V 1 . 
This is justified for a fully gapped Kondo insulator like 
YbBi2 where Vk 7 does not vanish along lines in k space. 
Furthermore to use only a minimum set of model param- 
eters, we replace V 1 = V T , m by V r = (l/2)(E m l^r, m | 2 )^ 
which is the average over each set of quartet states. We 
use a nearest-neighbor tight binding model with hopping 
t for the conduction electron bands e^. The spectral 
function of the experimental dynamical susceptibility of 
YbBi2 exhibits two sharp peaksQ- Therefore it is essen- 
tial that the two CEF quartets have two different average 
hybridization energies Vp (T = 1,2). 

The mean-field approximation to Eq. |T]) is obtained 
by taking b = (bi). Minimizing the ground state energy 
with respect to b and the Lagrange multiplier A(, leads to 
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where the following expectation values are introduced 
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Ekl/L^), ™ 7 = ^-E k ( d L d k 7 ) and 



J = 



Ek(/k 7 /ki)- 111 Eq. ©, n is the density of holes 
per site which defines the chemical potential /i. The 
mean-field Hamiltonian can be diagonalized. One ob- 
tains, H M f = Ek 7 .± ^ 7 ,±( k ) a k 7 ,± a k 7 ,±> wner e the hy- 



bridized bands have energies -E 7j ±(k) = 



± 



4? 7 2 



which are still fourfold (m=l-4) de- 



generate. The corresponding 4f- weight functions of these 
quasiparticle bands are given by A^ ± (k) — A^ T (k) = 

also obtain 



-e k ) 2 +4K? 



From the mean-field solution we 
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£ T;+ (k) - £ 7 ,_(k) : 
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?g(fc)/(B 7 ,±(k)), 



(3) 



k,± 



where = 



Af,. In the zero 



and e/ T 

temperature limit, T = 0, the upper bands are empty. 
Then the Fermi functions reduce to /(-E 7i +(k)) = 0, and 
/(-E 7 ,_(k)) = ifij (S-vVt = Under the condition, 

n = 2 or y> 7 = 1, which holds as long as the chemical 
potential is within the hybridization gap, we obtain the 
following mean-field equations from Eqs. ([2]): 



_ Vj? V(^- f 7 r ) 2 +4F r 2 + ^-e7 r 



r=i,2 



(D + e fT y+AV 2 -D-i fT 



r=i,2 

Here e/ 2 = e fl + A 2 , U 2 = Vi + 5V = 6(Vi + <5F), and 
D = 6t is half the conduction band width. The den- 
sity of states of the conduction band is assumed to be 
rectangular (g(e) — 1/2D; —D < e < D and zero 
otherwise). By solving the set of equations numerically 
one can find the tf t and b values. In order to be in 
the Kondo limit and have an insulating state with small 
hybridization gap the parameters should fulfill the con- 
dition A 2 < SV < Vi, | e/ |< D. In the absence of CEF 
effects, by choosing e/ = — 0.75i, V\ = t, SV = 0, A 2 = 
we found ef 1 = 0.05i and b — 0.30 from the mean-field 
solutions which will be used in Fig. [TJ 

The dynamic magnetic susceptibility is calculated 
within RPA approximation. Since we have two CEF 
quartets the spin response has the matrix form x(q, w ) = 
[/ — J q xo(ci, _1 Xo(cii w )i where the unit matrix I, the 




FIG. 3: The imaginary part of the susceptibility for the nonin- 
teracting ( Jq = 0) and interacting case ( Jq 1 = Jq 2 = 0.065t) 
for V% = t, q = Q and SV = 0.13t. The Jq^ are slightly sub- 
critical leading to a finite intrinsic resonance line width. 



interaction J q and the non-interacting quasiparticle sus- 
ceptibility Xo(q> <*0 are 2x2 matrices in the CEF quartet 
index r = 1,2. 

The exchange interaction J q between quasiparticles 
is assumed to be peaked at the AF wave vector Q = 
(ir, 7r,7r), i.e., the L-point because there the most pro- 
nounced magnetic response is found. In principle J q may 
be calculated to order (1/N 2 ) [liED] but this is strongly 
model dependent. We choose to parameterize J q in a sim- 
ple way: The interaction is peaked at Q or Q - 

r 

it has the Lorcntzian form J„ 



'q rr , 



(q-Q/+r 



\/3n and 

J Qrr> ' 



where Tq has the meaning of an inverse AF correlation 
length. Each element of the irreducible susceptibility ma- 



trix is calculated from the quasiparticle states as 21|: 



*r(q.") = £4 ± (k + q )4, T (k) 



k.± 



/(£r,±(k + q))-/(£r-. T (k)) 
£r >T (k)-£r',±(k + q)-u; 



(4) 



The non-diagonal elements of the interaction matrix cor- 
responding to interactions of quasiparticles with differ- 
ent CEF symmetry are neglected, implying Jq rr , = 
Jr(q)<5rr'- Then the RPA susceptibility is simply a sum 
of two contributions x rr (q> w ) from the two sets of hy- 
bridized bands: 



X(q,w)=£[l- J r (q) X £ r (q, w )]- 



L ^ r (q,^). (5) 



We now discuss the results of numerical calculations 
based on the previous analysis. In Fig. ([I]) we have plot- 
ted the real and imaginary part of \o T (q> w ) without 
CEF splitting (A 2 = SV = 0) versus energy for wave 
vectors q = and q = Q. One notices that Imxg r (Q, uj) 
has a strong low-energy peak due to a small indirect gap 
while Imxg 7 (0,o;) has a small peak at much higher en- 
ergy due to a large direct gap. The broad structure of 
the former is due to noninteracting single-particle exci- 
tations and the q and uj dependence is depicted in the 
inset of Fig. [3J 
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FIG. 4: Contour plot of imaginary part of RPA dynamical 
susceptibility with Lorentzian interaction Jr (q) for CEF-split 
quasiparticle bands, ef ± — 0.08i and b — 0.41. Here V\ = 
t and V 2 = Vi + 0.13t (J Ql = 0.125t, Jq 2 = 0.143i and 
Yq = 2 which satisfy the resonance condition in RPA formula: 
Jr(Q) = l/i?e[xo r (Q,^r)]), in direction (1,1,1). The peaks 
at the zone boundary (q = Q) appear at u\ — 0.047t and 
Lo 2 = 0.063/. By choosing t = 0.32eV (D = 1.92 eV) then 
uji — 15meV and L02 = 20meV which are comparable with 
experimental results. The inset shows the contour plot of 
imaginary part of dynamical susceptibility of noninteracting 
degenerate bands for comparison (Vi = It, SV = 0.13t, V\ = 
0.41t, and Jq = 0) in the direction (1,1,1). The color scale of 
the inset is the same as Fig[2] 



The density of states for noninteracting quasiparti- 
cles p 7 (w) = -J- J2u ± — ^7,±(k)) including the CEF 
splitting for the two sets of bands with A2 = O.Oli and 
SV = 0.13/ is plotted in the inset of Fig. (TJJ. The two 
hybridization gaps are different due to a finite SV. How- 
ever the latter is kept small enough to ensure that the 
chemical potential is within the gap. 

When the AF interaction Jr(q) is turned on, the 
imaginary part becomes Imx rr (q, uj) = F(ap, ?7r)/Jr(q) 
where a r = J r (q)Imxo r (q, w), r) F = 1 - 
Jr(q)-Rexo r (q, w), and F(cv r ,??r) = "r/(?7r + a r)- In 
that case the spectrum for q = Q moves to lower en- 
ergies and a narrow double-peak structure, i.e., the col- 
lective spin resonance excitations appear. Their energies 
are determined by the solution of rjp — 0. If they 
are lying within the indirect hybridization gap one has 
ar - * 0. Then the spectral function is a delta function 
nS(r]r)/ 'Jr(q) at the resonance energy ui r . The disper- 
sion of the resonance, is determined by the real part of 
Xo r (q, presented in Fig. [21 The plot shows that for 
q < Q the maximum of the spectral function follows a 
ridge which decreases in height and bends to higher en- 
ergy. This turns into an upward dispersion of the res- 
onance pole. Its endpoint in the BZ is limited by the 
extension of the ridge in Re\o- The latter is fixed for 
the simple hybridization band model used here. A more 



realistic band model might give a larger extension than 
the one seen in Fig. [2] 

Due to the CEF effect the f-levels split into two 
(pseudo-) quartets (A2 > 0) which hybridize differently. 
For SV > the resonance w^ =2 associated with the 7 = 2 
hybridized bands moves to higher energy and a second 
peak in addition to the one at appears in the spec- 

tral function. This is clearly seen in Fig. [3] where the 
CEF split resonance peaks at q = Q appear around the 
threshold energy of the non- interacting continuum states. 
In this figure we use subcritical values for the interaction 
constants. Therefore the resonance peaks are right above 
the continuum threshold and have a finite intrinsic line 
width. 

If the interaction constants are slightly increased the 
resonances move below the continuum and turn into true 
spin exciton poles without intrinsic line widths (within 
RPA). Their dispersion is shown in the main panel of 
Fig. [U Away from the L-point (Q = \/3tt) they dis- 
perse upwards and merge into the continuum. We iden- 
tify these spin resonance modes with the observed ex- 
perimental peaks at 15 and 20 meV 0, @] and we have 
chosen parameters such that their energy splitting and 
dispersion are reproduced. Our numerical calculations 
show that the best fit to experiments is obtained for 
SV = 0.13* where e fl = 0.08i and b = 0.41. These spin 
exciton peaks separate with increasing CEF splitting A2 
and hybridization energy difference SV. Therefore the 
influence of the latter is strong since it directly affects 
the hybridization gap and hence the noninteracting sus- 
ceptibility and resonance condition. We note that an in- 
crease in Jq (or a decrease of the hybridsation gap) will 
lead to a decrease of the spin exciton mode frequencies 
at Q. For Jq x = 0.179/; the lowest mode would become 
soft. This softening signifies the instability of the para- 
magnetic state and the onset of AF order in a Kondo 
semiconductor. This is not observed in YbBi2 at am- 
bient pressure. We suggest that an investigation of the 
pressure dependence of spin exciton mode frequencies at 
Q would give important clues how close YbBi2 is to 
AF order. Finally we mention that our present model 
does not include the broad excitations at 38 meV. As 
has been suggested in Ref. [l3] it might be due to a con- 
tinuum of additional band states which do not take place 
in the resonance formation. Their inclusion would re- 
quire a multi-orbital conduction band model. 
We thank P. A. Alekseev and I. Eremin for helpful dis- 
cussions. 
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